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Abstract: 

We present a study of the finite-volume two-pion matrix elements and correlation functions of the / = 
scalar operator, in full and partially quenched QCD, at one-loop order in chiral perturbation theory. In 
partially quenched QCD, when the sea and valence light quark masses are not equal, the lack of unitarity 
leads to the same inconsistencies as in quenched QCD and the matrix elements cannot be determined. 
It is possible, however, to overcome this problem by requiring the masses of the valence and sea quarks 
to be equal for the u and d quarks while keeping the strange quark (s) quenched (or partially quenched), 
but only in the kinematic region where the two-pion energy is below the two-kaon threshold. Although 
our results are obtained at NLO in chiral perturbation theory, they are more general and are also valid 
for non-leptonic kaon decays (we also study the matrix elements of (8, 1) operators, such as the QCD 
penguin operator Qe)- We point out that even in full QCD, where any problems caused by the lack of 
unitarity are clearly absent, there are practical difficulties in general, caused by the fact that finite- volume 
energy eigenstates are linear combination of two-pion, two-kaon and two-rj states. Our work implies that 
extracting AI = 1/2, K —> tttt decay amplitudes from simulations with ms = rnd,u is not possible in 
partially quenched QCD (and is very difhcult in full QCD). 



PACS numbers: 11.15.Ha,12.38.Gc,12.15Ff 



1 Introduction 



Several methods have been proposed to compute non-leptonic kaon decay amphtudes in lattice QCD. In 
principle one could determine physical amplitudes, including the final state interaction (FSI) phase shifts, 
from matrix elements of the relevant operators of the effective Hamiltonian, computed on the lattice in 
full QCD, with realistic quark masses, in a finite volume [1, 2]. This programme requires, however, 
computer resources which are not available at present, and which will not be available in the foreseeable 
future. For this reason several strategies for the determination of the matrix elements at next-to-leading 
order (NLO) in the chiral expansion have been presented in the literature [3]- [6]. The general idea is 
to compute matrix elements at unphysical quark masses and for a range of energies and momenta for 
the mesons in order to determine all the low-energy constants required at NLO in Chiral Perturbation 
Theory (xPT). xPT at NLO can then be used to obtain the physical matrix elements. 

A number of convenient choices have been suggested for the unphysical kinematics at which the matrix 
elements can be computed directly, and from which the complete set of couplings (low-energy constants) 
of the weak chiral Lagrangian at NLO in xPT can subsequently be determined. The original proposal 
was to compute the K tttt matrix elements with the kaon and one of the final state pions at rest, and 
to vary the momentum of the second pion (SPQR kinematics) [3]. Laiho and Soni (LS) on the other 
hand, propose to combine the K ^ 0, K ^ tt and — > tttt amplitudes, with the two pions at rest [5, 6]. 
In this case one has to isolate only the contribution from the two-pion state with the smallest energy in 
the correlation functions. LS also argue that this method can be used in partially quenched two-flavour 
QCD (PQ2), provided that the mass of the "sea" quarks in the loops, ms, is the same as the "valence" 
light quark mass, my = mu^d (the strange quark is quenched) [6]. 

In this paper we present a study of two-pion amplitudes in the / = channel and the corresponding 
correlation functions in a finite volume in partially quenched QCD. The results discussed below are 
important for recent attempts to develop new strategies for computing K tttt decay amplitudes using 
numerical simulations on the lattice. We also comment on subtleties in the evaluation of two-pion matrix 
elements in full QCD. 

We have previously used xPT at NLO to study the determination of matrix elements with two-pion 
external states in quenched QCD both for / = 2 final states (for which the difficulties discussed below 
are not present) [4] and for J = final states [7]. For / = final states, we had found that the absence 
of unitarity in quenched QCD means that the calculation of the amplitudes is not possible [7]. The 
LS proposal has prompted us to extend these studies to partially quenched theories. In this paper we 
confirm that the inconsistencies of the quenched theory, and in particular those resulting from the mixing 
of hadronic and ghost states, are also generally present in partially quenched QCD, unless the masses of 
the sea and valence quarks are degenerate. In PQ2, the masses of the sea and valence u and d quarks are 
indeed degenerate, the rj' is heavy and decouples from the low energy Lagrangian, nevertheless we show 
that the lack of unitarity prevents the extraction of "physical" amplitudes from finite volume correlation 
functions when the two-pion energy Wjr > W™"" ~ 2Mk^ where W™^ is the smallest two- kaon energy in 
the finite volume, corresponding to the two-kaon threshold (throughout this paper we assume that is 
sufficiently small that inelasticity effects due to states with more than two mesons can be neglected). We 
illustrate this point with explicit calculations in xPT of the correlation function of the scalar operator S 
with two pion fields (we had used the same example in our study of quenched QCD [7]). In addition we 
also consider (8, 1) operators, such as the QCD penguin operator Qq whose matrix elements are important 
in the theoretical prediction for e' /e, and confirm the same features. 

Our work clarifies the conditions necessary to evaluate matrix elements with two-pion external states, 
at NLO in xPT, in partially quenched QCD. In particular, in PQ2 the necessary condition < VF^'" 
means that it is not possible to extract the K ^ tttt matrix elements when ms = my = mg, where rus 
is the mass of the strange quark. Since these matrix elements are needed for the implementation of the 
LS proposal [6], we conclude that it cannot be applied in PQ2. 

Our work also shows that, in spite of the difficulties mentioned above, the two-flavour partially 
quenched theory, PQ2, can still be useful at least at NLO in the chiral expansion. We show that a 
consistent extraction of the "physical" matrix elements is obtained if we work with rriK > {rUs > 
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my = nis), at an energy of the two-pion state in the finite volume, Wt^, which satisfies tlie Lellouch- 
Liischer (LL) elasticity condition W„ < WjJ**" ~ 2Mk- The consistency of PQ2 with this choice of the 
parameters opens new avenues for the calculation of the A/ =1/2 rule and s' / s amplitudes from K tttt 
matrix elements following the proposal of refs. [3, 4]. For these transitions it would be very difficult in 
practice to determine the amplitudes from a full QCD three-flavour computation of the K — > tttt matrix 
elements, with the mass of the strange quark different from the light quark masses. A two flavour 
unquenched calculation, although difficult, is nevertheless, feasible in the near future. With the required 
kinematical conditions however, one loses the practical advantage of using only pions at rest ^. Below, we 
will also comment on some of the subtleties in the extraction of the K ^ ww amplitude with tuk = 
in fuU QCD. 

We make the further important, if somewhat disappointing, observation. In ref. [8] it had been 
suggested that it is possible to obtain exact physical information about full QCD from three-flavour 
partially quenched QCD (PQ3), i.e. from unquenched simulations with three flavours at ms 7^ my- 
The practical advantage of PQ3 is that it is computationally less expensive to vary, and in particular 
to reduce, the masses of the valence quarks. Having determined the low-energy constants of xPT in 
PQ3, one can then extrapolate to the physical point, ms = my- We have explicitly checked that for 
the cases at hand with two-pion final states (the matrix element of the scalar operator or kaon decay 
amplitudes) this suggestion does not work and one has the same difficulties as in PQ2: in infinite volume 
the amplitude is singular at threshold, there is mixing of hadronic and ghost states and the finite- volume 
correlation functions are plagued by terms which grow linearly or cubically with the volume thus making 
the extraction of physical amplitudes impossible. Moreover, even at fixed finite volume, the presence of 
terms which depend quadratically or cubically on the time distances prevents the reliable determination 
of the matrix element. Perhaps for some simple quantities, such as the leptonic decay constants of mesons 
or semileptonic form factors, the strategy proposed in [8] can be used. From the above considerations, 
we conclude that this is not true in general. 

The plan for the remainder of this paper is as follows. In the following section we discuss the extraction 
of the matrix elements in full QCD. We show that although the theory is unitary, for the degenerate case 
m-s = rn,ii = m-u the fact that in a finite volume the energy eigcnstates are linear combinations of two-pion, 
two-kaon and two-77 states makes the extraction of two-pion matrix elements subtle and difficult. We 
then turn to partially quenched QCD for which the explicit results at NLO in xPT for both the scalar 
operator and for (8,1) operators are presented in the appendix. In section 3 we study the implications 
of these results for the determination of two-pion matrix elements in partially quenched QCD. Section 4 
contains our conclusions. 

2 Two-pion matrix elements in full QCD with uik > ttIt^ and 
rriK = rriT, 

In this section we consider full QCD and present the finite volume correlation function of the scalar 
density operator S = uu + dd + ss 



at NLO in the chiral expansion. We take Q < t2 < t\. iTqit) is the Fourier transform of an interpolating 
operator for the pion {'K{t,x)) 



We only discuss correlation functions and amplitudes in the centre-of-mass frame of the two final-state 
mesons. 

We focus only on those terms which concern the final state interactions (FSI). These are the terms 
which generate finite- volume corrections to the matrix element which decrease as powers of the volume and 

^Alternatively, one can work with rriK 7^ W-n and additional ultraviolet power divergences have to be subtracted in the 
matrix element. 



(0|<-(ti)7r^r(t2)^(0)|0) 



(1) 




(2) 
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which give the shift in the two-meson energy [9]. The remaining terms are only subject to exponentially 
suppressed finite-volume corrections and are not relevant in this discussion. The complete expressions 
for the correlation functions at one- loop order in xPT can be found in ref. [7]. 

At one-loop order in chiral perturbation theory, for an arbitrary value of the two-pion energy, keeping 
only those terms which are relevant for the energy shifts and finite-volume corrections to the matrix 
elements, we find: 

„-Bti p-Et2 / o \ 

(0|<^(ti)v(t2)5(0)|0) = (-J2) [l-.- + h{h,t2)] , (3) 

where E = \/(p' + and / is the pseudoscalar decay constant at lowest order in the chiral expansion. 
We write Ib{ti,t2) in the form 

Ib{ti,t2) = '^^Y. {d+{w^)P{w^) + c+{wk)P{wk) + ^^Pi^n)) + • • • > (4) 
where the volume is a cube of size L^, and 

1 _ g2(E-M))t2 
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^^'^^ " 2iE-w) ' w,^^Jk^+m-, 
c+{w) = =T-7^ , d+{w) = 2c+{w) - 



Evaluating the sum, we find 

-Eti ^ — Etg 



(0K^(ii)7rr(t2)5(0)|0) = (--pj [l + ReMoo)+T(t2)] , (5) 



where ^tx> is the infinite- volume one-loop correction to the amplitude [7] and 

^fe) - - - - - "'.)} + - . («) 

V = J2k-w=E ^^'^ terms which are proportional to t2 are exhibited here. 

The energy shift (AW) can readily be extracted from the terms proportional to t2- For E < ttik < nin 
only the first term on the right-hand side of eq. (6) contributes, giving 

AW = -(2-^]^ (7) 

V 4Ey pL^ ^ ' 

For energies such that E > rrirj^K = on the other hand we find 

The case E < uik < rrirj was discussed in detail in ref. [7] and we refer the reader to this reference for 
details. We only remark that in this case only P{w^) in eq.(4) contains terms with vanishing denominators 
(and hence terms proportional to 1,2) ■ Since we were considering energies in the region E < m.K < in 
ref. [7], the terms proportional to 9(E — ttik) and 9{E ^ m,,^) were omitted in equation (23) of that paper, 
but here we also want to study the case tuk — rn,, ~ < E, for which P{wk) and P(w^) also have 
vanishing denominators, so that the energy shift and the finite- volume corrections to the matrix element 
are modified. This result is in agreement with the Liischer quantization condition [9], because the FSI do 
depend on whether we are above or below the two-kaon (and two-Tj) threshold. The consistency can be 
readily checked by computing the one-loop expression of the relevant matrix element in infinite volume 
in Minkowski space: 



{7r-{q) TT+i-q) \ S \0) = -j^ [1 + A^] = -j^ 



(9) 
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where 



h = {ml - 2s)A{m„) - sA{mK) - -g^^lm^) + • • • • (10) 
In the above expressions, s = (p^- +p^+)^ is the square of the two-pion centre-of-mass energy and 



.m- /, / 1 + V^~^^ \ . . .m^' 



Aim) = \/ 1 - 4 log \ -inell- 4 . (11) 



s 



-sJT- 



When E < mx < m^j, the s-wave phase shift, which is obtained from the coefficient of the imaginary 
part of lb, is given by the first term of the r.h.s. of eq. (10) ^: 



2s — ml 
167r/2 V ^ ^ s 



^(«) = ^7r72^\/l-4— • (12) 



When E > mx = mj^ = m^ on the other hand, all the terms in eq. (10) contribute and we obtain 



From the result for the phase shift in cq. (13), wc can extract the energy shift by using the Liischer 
quantization formula [9]. The energy shift obtained in this way agrees of course with the result in eq. (8), 
obtained from a calculation of the correlation function in the finite volume. 

For the degenerate case E > = m,r = m,,, the energy shifts discussed here correspond to 

finite- volume eigenstates which arc two- meson SU(i) singlet states [S is an SU{i) singlet operator). A 
convenient procedure for the evaluation of the matrix element is to compute the correlation function of 
S with the SU(3) singlet two- meson operator {ipip)^, 

S-M, ^ (0|(^^_,-(<i)^f(i2))i5(0)|0), (14) 

where (for compactness of notation, in the following we suppress the labels t\^2 and q) 

, . I 11/3" 

^^^''i " 7!^^^ 272''"^^ 2V 2''''' ^^^^ 

and 

KK = ^ {K°K° + K+K-) , 'JT'JT=^ (7r°7r° + ^2^+^-) . (16) 
v2 v3 ^ ' 

{(p<f>)^ creates or annihilates the SU(3) singlet state |1). which is a combination of two-pion, two-kaon 
and two-?7 states with the same flavour structure as the corresponding operator 

IA-A-) = -i= + , M = ;^(|A») + ^|,+^-)) . (17) 

All the meson states are s-wave, 7 = states. Similar expressions can be written for the octet and 27-plet 
operators considered below. Computing in addition the {tpifi) ^ - {tfifi) ^ correlation function 

(0|(<^L,^(ti)^,-(i2))i(4(-t2)<^-«(-ii))i|0) , (18) 



^For s > 4ml, we have Arg [(7r-(g)7r+(-q)|5|0)] = S{s). 
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one obtains the finite volume matrix element |(1|5|0)|fv by dividing the correlation function (14) by the 
square root of the four-point function (18) at large time distances, following the procedure explained in 
ref. [7] 

l-ITH , '°I(«'U('.W..)).^(0)|0) 

V(0|(v^,-((l)Vj(<2)),(»'J(-t2)«)-f(-<l)),|0) 

The TTTT matri x ele ment is then simply obtained by using the appropriate Clebsh-Gordan coefficient, 

(ttttIS'IO) = 2y273(l|5'|0). 

For the sake of presentation, it was convenient to discuss the extraction of the matrix element using 
the ratio on the right-hand side of eq. (19), with the correlation function S- (w) i in the numerator. 
Note however, that since S is an 5J7(3)-singlet, we could equally well have used the correlation function 
S - (tttt) . 

In the case of (8, 1) or (8, 8) operators the simple procedure outlined above does not work because 
these operators have non- vanishing matrix elements between the kaon state, \K^), and three different 
S-wave two- meson states with 1 = and 1^=0. The three states are the state |1) defined above and the 
octet and 27-plet states (created by the corresponding {w)^ and {ifip)^^ operators) 

[97 1 

In a finite volume the three two-meson states (|1), |8) and |27)) acquire different energies, and thus they 
appear in the correlation function with different exponentials in time. Thus, in order to obtain the matrix 
element {ttttIQqIK'^), for example, one has to disentangle the different contributions, and extract all the 
matrix elements (1, 8, 27\Qe\K°) and (1, 8, 27|1, 8, 27) by studying the K°-Qe-{(p(p) . and the {(fiip) ^-{^pLp) . 
correlators. Then a suitable combination of (1, 8, 27|Q6|A'") will give the required (TTTrlQel^") ampUtude. 
A further complication is that all the matrix elements should be computed at the same two-meson energy, 
which must be kept fixed in the infinite volume limit [1, 2]. Although possible in principle, the procedure 
appears to be very complicated to implement in practice. 




3 Two-pion matrix elements in partially quenched QCD 

In this section we discuss the evaluation of the same matrix elements, but now in partially quenched 
QCD. We consider the two cases: QCD with three sea-quark flavours for generic values of ms, my 
and nis (PQ3), and QCD with two light sea-quark flavours also with generic combinations of the quark 
masses (PQ2). These calculations correspond to theories with SU{6\3) and SU{5\3) graded Lie groups 
respectively. We have also studied the S'J7(4|2) case, reaching the same physical conclusions, and for this 
reason this case will not be discussed explicitly. 

PQ2 is partially quenched two-flavour QCD, with a mass ms for the sea quarks, a mass my for the 
valence light quarks and a mass for the strange quark, which is always quenched. The rather lengthy 
formulae for ms ^ my are presented in the appendix for both the scalar operator and (8,1) operators. We 
explicitly conflrm the expectations [6] that when the sea-quark mass is different from the valence-quark 
mass the problems are the same as in the quenched case [7]. The double poles induce terms which seem 
to grow linearly and cubically with the volume as well as ones which are quadratic or cubic in the time 
distance, see eq. (27). The extraction of the physical amplitude is therefore not possible, at least with 
our current level of understanding. In this respect we flnd that it is euphemistic to state that "the final 
state corrections can be significant" [6] ; we simply do not know how to extract the infinite volume matrix 
element from the finite volume correlation function and the problem cannot be overcome by going to 
larger lattices. This pathology is not peculiar to PQ2 but is also present with the same symptoms in 
PQ3. 
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For ms = my, for a generic energy E and a generic strange quark mass nis we obtain in PQ2: 

-ELi -Et2 / o \ 

(0|<-(iiK^(t2)S^''(0)|0) = (-j^) [l... + ir{ti,t2)] , (21) 

where we choose = uu + dd + ss and the sum is only over the valence quarks, 

I^iM) = -^^^ (rf+K)PK) + Ic+{wk)P{wk) + ^^^^(«-)) + • • • (22) 

and rrigg = 2m|- — m^. When niK > {nis > my = ms) and taking the two-pion energy such that 
Wtt < W^™ (at this order in the chiral expansion this means 2E < 2mK), only the term proportional to 
P{wtt) generates the energy shift and the finite volume power corrections, which are exactly the same as 
in the unquenched two-flavour theory. In this respect we agree with the conclusions of LS for ttik = 2m„. 
More generally, one has finite volume corrections under control and the validity of the LL formula, as 
in two-flavour full QCD, whenever Wt^ < 2m k- Indeed, Wjr < 2m k is precisely the elasticity condition 
under which the LL formula applies. If wc violate this condition however, then major difliculties arise. 
Consider for example the case mx = = mj^: each of the three terms contributes to the energy-shift, 
and the flnal coefficient does not correspond to the correct result either in SU(2) or in SU(3). The reason 
for this is that the "quenched" states containing a strange quark now also contribute to the singular terms 
in the sum. The matrix elements now correspond to final states with a ghost component and therefore 
wc arc unable to extract the required two-pion amplitudes. In full QCD in sec. 2, we were able to use the 
SU(3) symmetry in the m„ = md = ms case to determine the TTTr-component in an SU(3)-singlet state. 
Here we are unable to carry out the analogous procedure. We note that one could avoid this diSiculty 
by using a scalar operator which projects only onto two-pion states, such as uu + dd. In this case the 
problem of the mixing with intermediate states containing strange quarks does not arise. This feature, 
however, does not apply to the (8, 1) or (8,8) operators of the weak effective Hamiltonian. 

The problems above, which are related to the absence of unitarity in partially quenched theories, are 
also present in PQ3 when W^r is above the threshold for any multi-hadron state containing "partially- 
quenched" quarks. Therefore in order to extract the AI = 1/2 K ^ tttt matrix elements with m„ = 
md = ms in PQ3, one is forced to work with the full QCD limit of the theory. 

4 Conclusions 

In this paper wc study the extraction of matrix elements with two-pion final states from simulations 
in finite volumes in full and partially quenched QCD. In full QCD, not surprisingly, we confirm that 
the matrix elements can be determined in principle, up to exponentially small finite-volume corrections. 
We point out however, that even in this case, for three degenerate quark flavours there arc considerable 
practical difficulties because the finite-volume energy eigenstates correspond to linear combinations of 
TTTT, KK and rjr] states. The determination of a matrix element with a two-pion final state with a fixed 
energy requires therefore the extraction of several finite-volume matrix elements at this energy. Each of 
these matrix elements will, in general, need to be computed from a simulation on a different volume. 

For PQ2 we find that it is not possible to determine K tttt matrix elements in general, at least 
with our present level of understanding. For two-pion energies below the two-kaon threshold however, 
the theory is a consistent unquenched two-fiavour one, the LL formula holds and it is possible to extract 
infinite volume physical amplitudes. Of course the value of the couplings of the effective weak chiral 
Lagrangian and the coeflacients of the logarithms which appear at one loop are not those of full QCD. 
One may argue, however, that the contribution to these couplings due to the sea-strange quark loops 
are suppressed by the massive kaon propagators when W^r << 2m k, since there are no singularities 
associated with the strange meson loops. If this is true, PQ2 can provide very useful physical information, 
although we are currently unable to estimate the systematic uncertainty. Given the apparent importance 
of unitarity in enabling the determination of physically meaningful results, we believe that it may be 
theoretically better to work with the unquenched two flavour theory rather than to weight the gauge 
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field configurations with fractional powers of the fermion determinant in order to mimic the three flavour 
theory. A study of this issue in chiral perturbation theory may be very useful to clarify the situation. 

The disappointing news is the requirement that the sea and valence quark masses should be equal. 
It had been considered very useful to work with fixed sea-quark masses and to vary (and in particular 
to reduce) the valence ones, in order to calibrate the lattice spacing and to extrapolate the results to 
the physical point. In this respect, our results cast serious doubts on the applicability of PQ3 to the 
determination of matrix elements with two-meson (or in general, multi-hadron) external states. PQ3 may 
however, prove useful for the determination of some other quantities. 

For PQ2 with the two-pion energy at or above the two-kaon threshold (more precisely for rriK > TUn 
{iris > = ms), with Wji > W^"^) with our present level of imdcrstanding it is not possible to 
extract the two-pion matrix elements. The difficulties are those already discussed for quenched QCD [7]. 
Unfortunately this implies that the proposal suggested by LS of combining K ^ K ^ -k and K ^ ttw 
amplitudes, with the two pions at rest, to determine the complete set of couplings of the weak chiral 
Lagrangian fails in PQ2. An ingredient in the LS proposal is the extraction of the two-pion matrix element 
with rn^ = niK, which violates the consistency condition given above. Of course the LS proposal may 
work in full QCD, but requires care in the extraction of the K mr matrix element when mx = m^. 
In this case the full two meson — + two meson correlation matrix must be studied, as explained in sec. 2. 
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A Finite- volume correlation functions in partially quenched QCD 

In this appendix we present the results for the correlation functions used in the extraction of two-pion 
matrix elements in partially quenched QCD, PQ2. The Euclidean partially quenched chiral Lagrangian 
in the supersymmetric formulation is [10] 

= y str [(a^SP^t) (d^J^Pi)] - :^str [E^^t^ + ^t^f^] _ ^^c&g -f a (9^$o) {d,,^,) , (23) 

where S**' is the graded extension of the standard non-linear Goldstone field E in the full chiral per- 
turbation theory, $0 is the super- 77' field and x — ^BqM, where A4 is the mass matrix and Bq = 
-{0\uu + dd\0)/f. 

The results reported here are obtained with the super-jy' integrated out (the limit mo — > 00), and are 
presented in terms of the following quantities: 

c±{w) = -r-^ ; cq{w) = -—^; rf±,o(w) = 2c±,o(w) - 



Ws = yk^+m'^g, myg^ Bn{mu + ms) and (24) 



Wi 



''vs ' '"ys 

mf - ^vs 
mj - m?- 



= ^Jp + m?^ (for i ^ S), 5s = {ml - m^^) , Bi. 



The finite volume scalar correlation function can be written in the following form {E = \/q'^ + ml): 

-Eti -Et2 / o \ 

{0\n+^^{h)nZ{t2)SP0{0)\0) = [-j^j [1 + /f + /r + /^(ti, ^2)] , (25) 
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where /f' and 7^' are the contributions from the wave-function renormalization and tadpole diagram 
respectively. Thoy give the infinite volume result up to exponentially suppressed corrections with the 
volume and are therefore not reported explicitly here. 

iritiM) = -2^2 i^ss + Akk + + Aq^ + Am] , (26) 

where 

A ^ 1 [ C+{ws) C-{ws) 2(E-ws)t-,. ( C+{ws) Co{ws) \ 

i3 4^ 1 2{E - ws) 2{E + ws) \2{E-ws) 2ws J 

k 

, 2Et2-2wsti ( Co(wg) _ C_(wg) 

2ws 2{ws + E) 



^ = _ X 1 / '^+i'^K) _ C-{wk) _ 2{E-WK)t2 ( C+{wk) Co{wk) 

2 \2{E - wk) 2{E + wk) \2{E - wk) 2wk 

k 



r,3 X]' 



T,3 m 



k 



2wK 2{wK + E)^ 



^\ 2{E - w^) 2{E + w^) 



_^2{E-w^)t2 ( d+{w^) - c+{w^) ^ do{w^) - co(w^) 



+e 



2{E - w^) 2w^ 

2Et2-2w^U ( - Co(w^) d-{w^) - C-(W^) 



2w^ 2{w^ + E) 

,2 10:2 I „,4 d2 



1 ^ ml{25l+wlBl^) f _^ l__ _ ( 1 , ^\ .27^ 

^L^L. 4£2y,6 \2{E-w^) 2{E + w^) \2{E - w„) ^ 2w„ J ^ ' 



_|_g2Et2-2w„ti 



1 



2w,: 2{w.„+E) 

1 ^ _ m|5|_ r 1 1 _ ^2iE-n,^)t2 [ I - [2{E - w^)]h _ 1 + [2w^]t2 

E^wl 1 [2{E - w^)]^ ^ [2{E + w^)]^ \ [2{E - w^)]^ [2w^]^ 



^2Et2-2w^t^ [ l + [2w^]tl _ 1 + [2{W^+E)]ti 

[2w^]^ [2K+^)]2 

1 



j_^m|^ r 1_ 

L3 E^wi \[2{E- 



)]3 [2{E + w^)r 



_,2(E-^^)t2 ( 1 - - w^)]t2 + [2{E- w^)]Hl/2 1 + [2w.]t2 + [2w^]Hll2 



[2{E~w^)Y ^ [2w.f 

,^2Et2-2^,t, ( \+\2w^\ti + \2w^ft\l2 _ 1 + [2K + + [2K + -E)]'^i/2 
V [2K + i;)]3 

= — V "^^^^^^^^ / ^ 1 g(2iS-«;^-«,j,)t2 1 ^ 1 

i3 ^ 2E'^Wt,Wss \2E-Wt,- Wss 2E + WTr+ Wss \2E-W.^- Wss Wtt + Wss , 

k 



_|_g2Et2-(«;,+t«sa)ti 



1 



Wtt+Wss Wss -\- Wt^ 2E ^ 

1 1 



iS Z-/ AE'^wl { 2{E - Wss) 2{E + Wss) \ 2{E - Wss) 2ws 

k 



j^^2Et2-2wssti 



1 1 



2wss 2{wss + E)^ 

Note that the presence of flavour-singlet double poles gives rise to terms proportional to powers of ti and 



8 



t2 and that these terms disappear when mg = my (when ds and B^k vanish). 

The same structure is present in the K ^ ww correlation functions of the (8,1) operator. We denote 
this operator by Q^^ (Qq is an example of a QCD operator which transforms as the (8,1) representation), 
and in the chiral effective theory it is chosen to be 



QP^I^ = str [Af a^E«a^Ef«t] , where 

Ae 




Af = ( " ) • (28) 



The correlation function is given by 



(0|<~(ti)^^(t2)Q^|,,)(0)if2(t^)|0) = 



e-^*i g-_Ei2 ^rriKtK / gj 

~~2E 2E 2mK \~P 



T{E,mK,m„) - — 
2/^ 




where the ellipses represent non-singular contributions, the sum runs over all mesons and 

,^ , 2E^ + Emi-ml 
t{E, mi,m2) = . 

To this order the relevant contributions to Aj can be written as 



^ E E ( 1 - e^-— E(2i. - - «^.)"^J + ... , (30) 



L^^^^{2E-Wi-wjy 



where again only the terms containing poles have been kept. The only non-zero coefficients are the 
following: 

C^KK = -T{WK,mK,mK)c+{WK) , 

^si = r{ws,mK,mvs)c+iws) , 



ml 
E^wi ' 



^(1) D D ml (Wss+W-K \ mlniKSsBKn 

CkJ = -B-^kBktv-p^ T ,mK,mK ] 



^(2) _ mlmK5sBKn{wss - w„) 



2E^wl ■ 

where the quantities on the right-hand side are defined in eq. (24). 
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